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Abstract: We investigate the gravitational lensing scenario due to Schwarzschild-like black 
hole surrounded by quintessence (Kiselev black hole). We work for the special case of Kise¬ 
lev black hole where we take the state parameter w q = — For the detailed derivation 
and analysis of the bending angle involved in the deflection of light, we discuss three special 
cases of Kiselev black hole: nonextreme, extreme and naked singularity. We also calculate 
the approximate bending angle and compare it with the exact bending angle. We found 
the relation of bending angles in the decreasing order as: naked singularity, extreme Kiselev 
black hole, nonextreme Kiselev black hole and Schwarzschild black hole. In the weak field 
approximation, we compute the position and total magnification of relativistic images as 
well. 

Keywords: Black hole; gravitational lensing; null-geodesics; quintessence; relativistic im¬ 
ages. 


I. INTRODUCTION 


Gravitational lensing (GL) signifies the deflection of electromagnetic waves. Light propagates 
in empty space along a straight line. The well-known theory of General Relativity (GR) predicts 
that light will be bent if an object with a certain gravitational field is interposed in the light path. 
In literature, GL has been used to study highly redshifted galaxies, quasars, supermassive black 
holes, exoplanets, dark matter candidates, primordial gravitational wave signatures, etc., jl] . In 
1801, Soldner was the first person who calculated the bending angle of light by using Newtonian 
Mechanics [2]. In 1911, Einstein derived the same Soldner’s result by using the equivalence prin¬ 
ciple and Minkowski metric, unaffected by gravity {3]. This marks the beginning of our modern 
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understanding of GL. In 1915, Einstein derived the new solar light deflection angle that was dou¬ 


ble from the previous value due to the effect of the spacetime curvature 


confirmed the prediction of Einstein dnnng tire solar eciipse Q. In 1937, Zwicfiy estimated the 
gravitational lens effect can be observed [6]. In 1979, Walsh, Weymann and Carswell used Zwicky’s 
work and discovered the first example of GL in which they obtained the first multiple images of a 
binary quasar (QSO 0957 + 561) jj]. 

In 1959, Darwin calculated the light deflection angle due to a strong gravitational field using 
the Schwarzschild metric 8|. Another significant work involved the deflection angle and intensities 
for the images formed due to the Schwarzschild black hole in terms of elliptic integrals of the first 


4|. Eddington in 1919, 


kind [9]. Considering the Schwarzschild black hole for the strong GL, Virbhadra and Ellis obtained 
the lens equation and introduced a method to calculate the bending angle. They also studied 
the lensing problem for the galactic supermassive black hole numerically jldf ]. While studying GL 
with the Schwarzschild black hole in the strong field limit, the bending angle was also evaluated 
analogous to the weak field limit. Besides the weak field limit of relativistic images, magnifications 


and critical curves formulas were also formulated [ll ]. Bozza treated the strong lensing phenomenon 
3 V a spherically symmetric black hole, where an infinite sequence of higher order images are formed 
uj] and later on extended for a spinning black hole [13]. One of the first important studies about a 
cosmological constant relativistic bending angle was done by Rindler and Ishak where they showed 
that for a Schwartzschild de Sitter geometry, the cosmological constant does not contribute to 
the bending angle [l4]. Another important application of relativistic bending angle techniques 
were used to determine a limit in the cosmological constant by using the bending of light through 
galaxies and clusters of galaxies jjj|. 

About two decades ago, a very important astronomical observation (using Super novae type la) 
suggested that the Universe is in a state of an accelerated expansion flfil. Tril]. This study was 
a revolution in physics and the dark energy was named to be responsible for this accelerating 
scenario. Cosmologists proposed different models in order to explain this strange behaviour of 
the Universe such as the ACMD model (with a state parameter of w = —1) or dynamic scalar 
fields 19]. The former uses the old idea of a cosmological constant introduced by Einstein 
several years ago but in a completely different way, 1 now interpreted like a reason to support the 


1 Einstein introduced a cosmological constant in his field equations to obtain a static universe. After some observa¬ 
tions that suggested that the Universe is expanding, Einstein thought that this constant was the worst mistake in 
his life. However, nowadays, this constant has been taken into account but using another physical interpretation 
related with dark energy 
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dark energy. However, this model has some problems like the so-called “cosmological constant 
problem” where the value of the cosmological constant differs about 10 120 orders of magnitude 


from the empirical value 


B- 


The second candidate for dark energy is a dynamic scalar field such 


as quintessence, phantoms, k-essence, etc. 


21-123]. Generally, a quintessence model has a state 


parameter w(t) = p(t)/p(t), where p(t) is the pressure and p(t) is the energy density that varies 
with time depending on the energy potential V(<h) and scalar field <f>. In addition, it is important 
to mention that the quintessence field is minimally coupled to gravity and the potential energy 
decreases as the field increases. This model is the simplest case without having theoretical problems 
like Laplacian instabilities or ghosts. For a more detailed review of the quintessence, see 124-1261] 
One important solution related to the quintessence model was discovered by Kiselev 27]]. The 
former solution physically describes a spherically symmetric and static exterior spacetime filled with 
a quintessence field, hence a nonvacuum solution. The Kiselev obtained the Schwarzschild-like and 
Reissner-Nordstrom-de Sitter BH’s solutions surrounded by the quintessence at the range of state 
parameter —1 < w q < — the Universe will accelerate with the quintessence, where w q is the ratio 
of pressure and energy density of quintessence. At w q = — 1, quintessence covers the cosmological 
constant A term and corresponds to the case of dark energy, while w q < — in a static coordinates 
quintessential state, reveals a de Sitter type outer horizon. In short, the solutions that corresponds 
to — 1 < w q < — i are asymptotically de Sitter. In this paper, we study the gravitational lensing 
due to a Kiselev black hole (KBH) where we choose the state parameter w q = — Due to this 
value, the solution will be a Schwarzschild-like (netural) black hole surrounded by quintessence 
[27]. In this paper, we considered three possibilities for KBH: two distinct horizons (nonextreme), 
unique horizon (extreme black hole) and no horizon (naked singularity). From the astrophysical 
point of view, it is a hard task to distinguish between the signatures and properties of black hole 


and naked singularities, however, GL can provide distinguishing signatures 


28]. 


The paper is structured as follows: In Sec. II, we study the geodesics and effective potential 
for nonextreme and naked singularity. In Sec. Ill, we discuss critical variables and equation of 
path for photons and calculate the relations between closest approach r Q and impact parameter b. 
In Sec. IV, we derive the bending angle in terms of elliptical integrals for both nonextreme KBH 
and naked singularity for different values of quintessence parameter a (discussed later) and then 
make a comparison with the bending angle for a Schwarzschild black hole. In Sec. V, we study 
the geodesics and effective potential for extreme KBH. In Sec. VI, we discuss critical variables 
and the equation of path for photons and calculate the relationship between the closest approach 
and impact parameter for the extreme lensing scenario. In Sec. VII, we calculate the bending 
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angle in terms of elliptical integrals for an extreme Kiselev black hole (EKBH) at a fixed value of 
a and compare it with the Schwarzschild bending angle as a reference. In Secs. VIII, IX, X, we 
use an alternative method for finding the bending angle to study the relativistic images. Finally 
we discuss our results in Sec. XI. We adopt the units c = G = 1. 


II. BASIC EQUATIONS FOR NULL GEODESICS IN KISELEV SPACETIME 


The equation of state parameter w„ for the quintessence scalar field <h is given by 


w„ = — = 


P q I$ 2 +V($)' 


Pg 


(1) 


where p q and p q are the pressure and energy density of the quintessence held defined in terms of 
the kinetic energy (^T 2 ) and potential energy 14(4*), respectively. Here, the overdot represents the 
differentiation with respect to cosmic time. 

Based on the above point of view, the geometry of a static sp 
surrounded by the quintessence (or Kiselev spacetime) is given by 


rerically symmetric black hole 


27J] 


1 


ds 2 = f(r)dt z — jr^dr 2, — r z d6 z — r z sin 2 ddcj) , 

fir) 


where 


f( S _ , 2 M a 

fv) 1 r r 3w q +l 


( 2 ) 


Here M is the mass of the black hole and u is the quintessence parameter (normalization factor) 
that is related to the energy density as follows J 27 ]: 

a 3 w n 


P q 2 7 - 3 ( 1 +™,) ‘ 

When w q approaches —1, the function /(r) for the metric (J^J) reduces to 

u \ 1 2fH 2 

fir) = 1- <rr 2 , 

r 


( 3 ) 


( 4 ) 


which is the Schwarzschild-de-Sitter black hole spacetime. For this case, the lensing phenomenon 
has been studied by Bakala and others [29143ll]. In this paper, our focus is on the special case 


w q = — which corresponds to the Schwarzschild-like black hole surrounded by quintessence. In 
this case the function /(r) becomes 

M. 


, . 2 M 

fir) = 1- err, 

r 


0 < a < 


8M J ’ 


( 5 ) 
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which can also be written as 


f( r ) 



r_)(r — r + ). 


( 6 ) 


The metric (J2J) becomes ill defined at r = 0, i.e., (goo —> oo) which gives a curvature singularity. 
For f(r) = 0, we get two fixed values of r, namely 


r + = 


1 + VI - 8Mcr 
2a 


r. 


1 - VI - 8 Ma 
2a 


(7) 


The region r = r_ corresponds to the black hole’s event horizon while r = r + represents the 
cosmological event horizon. Note that both r_ and r + are the two coordinate singularities in the 
metric ([2]). The coordinate singularities arise when 0 < a < However when a > both r+ 
and r_ become imaginary, giving a naked singularity. When a = 0, r_ becomes the Schwarzschild 
BH’s event horizon rfj = 2 M. 

The Lagrangian for a photon travelling in Kiselev spacetime is given by 


C = (l — llM- — ar ^J i 2 -_1- r 2 — r 2 0 2 — r 2 s i n 2 (8) 

r or 

Here dot represents the derivative with respect to A which is an affine parameter. We will work in 
an isotropic gravitational field, thus we can restrict the orbits of photons in the equatorial plane 
(0 = ■|). Hence, Eq. ([8]) becomes 

c = ( x - v- - ar ) p - 1 _A _ a / 2 - r ^ 2 - (9) 

By using the Euler-Lagrange equations for null geodesics, we get 


- _ dt 
d\ 



E 


1 - 
L 

3 ’ 


2M 

r 


— ar 


( 10 ) 

( 11 ) 


where E is the energy per unit mass and L is the angular momentum per unit mass. Using the 
null condition of the 4-velocity g^ y u^u v = 0 (where fj,, v = t, r, 9, 0) and lE = Vy- known as the 
4-velocity we get the equation of motion for photons, that is 


/1 

1 

/ 2 M \ 

where b = 

L 

A / TT7 —^ 

1 - ar). 


V b z 


\ r J 


E 


( 12 ) 


Here b is the impact parameter for photons of finite rest mass 


32], and it is the distance perpen¬ 


dicular from the centre of the black hole to the normal line on the ray of light intersecting the 
observer at infinity 33]. 
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The motion of geodesics is a force-free unaccelerated motion. In the presence of a gravitational 
field, photons experience gravitational force and this force comes due to the effective potential. 
Here, the effective potential for photons travelling in spacetime ([2]) is given by 


H eff = 



(13) 


Note that the effective potential has different values of a for nonextreme, extreme and naked 
singularity of KBH, i.e., for nonextreme 0 < a < g^j, for extreme a = while for naked 
singularity cr > . Here we discuss nonextreme and naked singularity cases and the extreme case 

will be discussed in Sec. V. When a = 0 then Eq. (fl3l) reduces to Schwarzschild BH’s effective 
potential, i.e., 



(14) 



FIG. 1: Effective potential V e s of photons as a function of distance r from black hole, setting M = 1. Top 
curve for Schwarzschild black hole, middle two curves for nonextreme while bottom two curves for naked 
singularity of KBH. 

In Fig. |H the effective potential 14s is plotted to study the behavior of photons near the 
considered spacetime ([2]) for different values of quintessence parameter a. We take M = 1 for 
plotting cr = | = 0.125 and the limits on a become for the nonextreme case 0 < a < 0.125, for the 
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extreme case a = 0.125 (discussed later in Sec. V), and for naked singularity a > 0.125. Hence 
(j = 0 corresponds to the Schwarzschild black hole, a = 0.06 and 0.1 corresponds to the nonextreme 
KBH. For these cases photons do not cross the horizon while at <7 = 0.14 and a = 0.15 photons 
cross the horizon. In each curve there is no minima. Therefore, there is no stable orbit for the 
photons, only an unstable orbit exists in each case which corresponds to the maximum value I4nax- 

III. CRITICAL VARIABLES AND THE EQUATION OF PATH FOR PHOTONS FOR 

KBH 


To find the radius of circular orbit of photons, we use the condition = 0 to obtain 

1 ± y/1 - 6 Mo 


r c ± = 


a 


(15) 


Here r c + is greater than the outer horizon r+ while r c _ lies between the inner and outer horizons 
(r_ < r c - < r_|_). The region of interest is between the horizons. Therefore, the radius of an 
unstable circular orbit for a photon is r c _ = r ps , also called the photon sphere. For the critical 
value of the photon sphere, conditions imposed on a are 0 < er < for the nonextreme and 
(7 > gjj for naked singularity. In the limit cr —>• 0 we get the radius of photon sphere r ps = 3 M for 
the Schwarzschild black hole. Now, we convert the equation of motion (1121) in terms of u = = We 
obtain the equation of path for photons 




(16) 


where 


B(u) = =r — u 2 (1 — 2 Mu — —). 
b z V u/ 


(17) 


For critical value of the closest approach, we put ^ = 0 [9j. Identifying this point of the closest 
approach as u = U 2 , from Eq. csd, we have 

1 


b 2 


= — 2Mu'2 — au2- 


(18) 


Substituting U 2 = — from Eq. (fl5l) in Eq. (fl8l) . we obtain the critical value of impact parameter 


for circular orbits 


bar — 


ps 


r ps - 2 M - oTp S ' 


(19) 
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The value of the impact parameter also imposes the same limits on the quintessence parameter 
<t, for both nonextreme and naked singularity of KBH as mentioned above. For a = 0, Eq. (1191) 
gives the impact parameter b^ c = for a Schwarzschild black hole. According to the circular 

orbit condition (setting B(u) = 0) and solving Eq. ([TTJ) . we get one real root u\ and two other 
roots U 2 and 7 / 3 , ( 7 t 3 > 112 > u\) which are 


r 0 — 2M — J(1 - 8 Ma)r'i + 4 Mr a - 12M 2 1 r a - 2M + - 8 Ma)r'j + 4Mr a - 12 M' 2 , N 

“l = - 1 -7T7-i U2 = —, U3 = -—-. (20) 


4 Mr 0 

Thus Eq. (1171) becomes 


4 Mr 0 


B{u ) = 2 M(u — u\){u — U 2 ){u — 113 ). 


( 21 ) 


Substituting Eq. (12T1) in (fl 6 t> yields 


du 


— = ±\/2M(u — u\)(u — U 2 ){u — U 3 ). 


( 22 ) 


In Eq. (25), the positive sign (+) shows that the angle (j>; changes more than 7 r, that is the 
photon trajectory is bent toward KBH and for the negative sign (—) the photon trajectory is bent 
away from KBH. For a ray of light, both r Q and b are obviously different from each other. Using 
Cardano’s method solving the cubic equation, 


r 6 a + crb 2 rl - b z r 0 + 2 Mb z = 0 , 


(23) 


the relation between b and r 0 is 


r Q = 2 


cr 2 6 4 + 3 b 2 


9 


■ cos 


3 C ° S 


-1 


2 (J 3 6 6 + 9ab 4 + 54 Mb 2 
6cj 2 6 4 + 186 2 


9 


a 2 b 4 + 3 b 2 J . 


ab 2 

"3 


- (24) 


At a = 0, it consistently reduces to the Schwarzschild black hole lensing case 


33| 


2 b 

r 0 = cos 
\/3 


1 _i /-3\/3 M\ 

L3 C ° S ^ C-6-J]' 


( 25 ) 
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FIG. 2: Closest approach r 0 as a function of impact parameter b {M = 1). We discuss here the relation 
between the closest approach r 0 and impact parameter b for KBH lensing cases: -nonextreme and naked- 
singularity- and compared it with a Schwarzschild black hole lensing case for different values of a. 

From Fig. [21 we observe that by increasing the value of 6, r a increases. In the region of the 
photon sphere a = [0, 0.1], r a depends on b from the quintessence parameter a. Moreover, as 
a increases, light moves closer to KBH and the closest approach r a decreases. Therefore, er = 0 
corresponds to a Schwarzschild black hole (taken as a reference) while a = 0.02 to er = 0.1 
correspond to the nonextreme KBH. Beyond the photon sphere (region where no horizon exists), 
i.e., a = 0.150, the light goes into the KBH, whereas r Q remains constant and naked singularity 
occurs. 


IV. BENDING ANGLE 


Suppose that a light ray comes from infinity (say — oo), reaches the black hole at r 0l and finally 
moves back to infinity (say + oo) that is the observer. Due to this change, the angular coordinate 
<f> is two times from infinity to r Q . The light ray deflects from a straight line path at the difference 
of 7r which results in the bending angle a fl 


a = 2 


I" r a d(j) 

Jo du 


du — 7T. 


(26) 
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If we substitute Eq. (1221) into Eq. (1261) , we obtain 


a = 2 


P 

Jo 


zdu — IT. 


(27) 


/0 a /2 M(u - Ui)(u - U2)(u - u 3 ) 

If we write Eq. (1271) in terms of complete elliptic integral 2 and an incomplete elliptic integral 3 we 
need to separate the integration limits into two parts: 


a = 


P 

J U\ 


1 


zdu — 


1 


zdu 


- 7r. (28) 


«l \/(ui - u)(u - u 2 )(u 3 - u) J U1 \J{u\ -u){u- U 2 )(li 3 ~ u) 

Here the integrals can be recognized in terms of a first kind of elliptical integral, where U3 > U2 > u\ 
351 ]. Hence 



M 


The integral variables can be defined as 


a = 2 /ir f ( f i- fc ) F(* 2 ,k) 


V u 3 ~ Ul V u 3 ~ Ui J 


— 7T. 


(29) 


'hi = 


7T _i r a — 2 M — a/( 1 — 8Mcr)r 2 + 4Mr 0 — 12M 2 

—, W 2 = sin \ - 

2 V r Q - 6 M - 

In the elliptical integral modulus k has a range 0 < |£:| 2 < 1, where 


(30) 


k = 


' 6 M — r Q + y/{ \ — 8 Mcr)r% + 4 Mr 0 — 12 M 2 


(31) 


2y /(1 - 8 Mcr)r 2 0 + AMr a - 12M 2 

Now F(^,k) = K{k ) defines a complete elliptical integral while F{^>,k) is an incomplete elliptic 
integral. By simplifying Eq. (1291) . an exact bending angle can be obtained: 


a = 4, 


y^l — 8 Ma)r% + 4Mr 0 — 12 M 2 


K{k) — F('h, k) 


— IT. 


(32) 


From the last expression, a can be deduced for nonextreme KBH under 0 < a < and for naked 
sin gula rity KBH under a > For a = 0, Eq. (|521) . reduces to the Schwarzschild bending angle 


a 


33]. 


2 The integral involving a rational function which contains square roots of cubic or quartic polynomials. Generally, 
here a definite cubic integrand that has a built-in command as 


de 


K (m ) = F {§ | m) = f * - 

y 1 —m sm z 9 

3 If 0 has the range — ^ < <p < 5 then F(<p 


*) = /„“ 


do 




m sin^ 9 
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FIG. 3: Bending angle is a function of impact parameter b. This is the case of nonextreme KBH lensing 
and its maximum deflection value depends on the quintessence parameter 0 < a < g (M = 1). Here the 
Schwarzschild case occurs at a = 0 while cr = 0.02 to 0.08 for nonextreme case. 

Figure [3j shows that the maximum deflection of light will occur at the critical value of the 
impact parameter b sc in Eq. (1171) . Below b sc there will be no deflection and above b sc , we will 
get a continuous deflection (light circulates around the black hole). Each single curve shows that 
by increasing the value of b, the bending angle decreases at different values of a. Nevertheless, 
originally when we increases the value of a, the critical value of the closest approach decreases 
since the light goes closer to the black hole. Similarly, the value of b (near the photon sphere 
where maximum deflection occurs) decreases and the bending angle increases. 

Figures [4] and [5] display the behavior of naked singularity. In Fig. 01 for any curve at short 
distances, as b increases the bending angle increases. In Fig. 01 for a long distance, as b increases 
the bending angle remains constant. However, when we observe the whole phenomena, we see that 
the bending angle also depends on a. As a increases, the bending angle decreases for both short 
and long ranges distances. Furthermore, when we compare the graph (Figs. 0]and[5]) of the naked 
singularity bending angle with the nonextreme and extreme bending angles graphs (Figs. [3] and 0}, 
we observe that naked singularity behaves opposite from nonextreme and extreme cases. 
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FIG. 4: Bending angle a as a function of b for naked singularity. At M = 1, a > I. 



FIG. 5: a as a function of b for a naked singularity. 

V. GRAVITATIONAL LENSING BY EXTREME KISELEV BLACK HOLE 

Extreme gravitational lensing is very amazing for some important phenomenona but it demands 
a great effort to be observed. In extreme gravitational lensing, where KBH is used as a lens, we 
need to discuss the bending of photons that pass very close to the lens and suffer a very large 
deflection. 
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For the extreme Kiselev black hole (EKBH) we have a = 1/8 M, thus the function f(r) becomes 

2 M 


fir) = 1 - 


r 

8 M' 


(33) 


This is an EKBH case for which /(r) = 0 gives rfj = 4 M known as a degenerate solution (single 
horizon). This value is twice the Schwarzschild black hole horizon, so it can be written as rj/ = 2rj|. 
Repeating the same procedure of Sec. II, for a = ^ we obtain the effective potential 


V e — — 

K eff- r 2 - 


L 2 2 ML 2 L 2 


8Mr’ 


(34) 


where the first term is related to the centrifugal potential. The second term represents the rel¬ 
ativistic correction due to general relativity. The third term arises due to the fact that EKBH 
geometry depends on a parameter a = g^y. Due to the effect of this potential, we can see the 
behavior of a photon surrounding by the EKBH. 



FIG. 6: Effective potential V^ s is shown as a function of distance r taking for extreme Kiselev lensing 
phenomenon. Observe that there is no minima (have no stable orbit) and only one maximum Enax, an 
unstable orbit that exists which corresponds to E/( ax . Schwarzschild’s effective potential is taken as a 
reference (cr = 0). 
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VI. EQUATION OF PATH AND CRITICAL VALUES FOR EKBH 

Substituting a = in Eq. (fT2l) . we obtain the first order nonlinear differential equation for 
path 

(S) 2 - B ' ( “ )=0 - (35) 

where 

B'W = 4-» 2 (l-2M«- s U). (36) 

In Eq. (1361) we need o apply the circular orbit condition. This condition gives a cubic equation 
that has one real root u\ < 0 and two distinct positive roots such that > 0. The roots are 


e _ rg - 2M - 2 y/(r% - 3M)M e _ 1 
U1 AMr e a ’ U2 r e 0 

Therefore, Eq. (|36l) can be rewritten as 


e _ r% - 2M + 2 ^{r% - 3 M)M 

iQ - - 

3 4 Mr e 0 


(37) 


B e (u ) = 2 M(u — u\)(u — u\){u — u%). 


(38) 


If we replace again this equation into the equation of path, Eq. f|3o[) . we obtain 

du 1 

d(j> y/2M(u — u\)(u — u|)(u — u§) 

In the limit u = 0 (r —>• oo), Eq. (1351) gives 


(39) 


0 

u = — + constant, 
b 


(40) 


For the critical value of the closest approach (radius of photon sphere r Q ), applying the second 
circular orbit condition ^\ u= j_ = 0, and then the condition dB d \ y | _i_ = 0 in Eq. f[35|) . we get 

= AM and r®_ = 12M. Here, r®, = r|j gives a degenerate solution (with 6 = 0) whereas 
r c- = r ps gives the photon sphere. Now, by putting the value of b^ c into Eq. (f35j) and using the 
condition of circular orbit B e (u) = 0, we get the critical value of the impact parameter, which is 
bg C = 6 y/GM. For EKBH the relation between r Q and 6 is 


r 


e 

o 


bVb 2 + 192 M 2 
12M 


cos 


.3 C ° S 


{ 


(6 4 + 2886 2 + 13824) ~ 
6 2 (6 2 + 192M 2 )§ J- 


b 2 

24 M' 


(41) 
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FIG. 7: Closest approach r 0 as a function of the impact parameter b for the EKBH. We see that by increasing 
the value of the impact parameter b the closest approach r 0 increases. Schwarzschild black hole case (a = 0) 
is taken as reference while for EKBH we take cr = 0.125 with M = 1. 

VII. BENDING ANGLE FOR EXTREME KISELEV BLACK HOLE 


The bending angle for the extreme Kiselev black hole (EKBH) can be obtained by putting Eq. 

(f39l) into (|26j) where r Q —>• r%. Doing this we obtain 

1 

1 


d e = 2 


f 

Jo 


zdu — 7r. 


(42) 


/ o y / 2M(rt — uf)(u — n|)(ri — u%) 

We can decompose the limits and convert the integral into complete and incomplete elliptical 
integral forms as follows 


1 


zdu — 


1 


VK -u)(u-u$)(ul -u) 


M L Ju\ yj{u\ - u){u - uQ(u% - u) 

Both integrals can be recognized in terms of first kind of elliptical integral 

e 
2 

2F{*lk' 


zdu 


— 7 r. 


(43) 


35y, where the integrand 


has the condition > u 3 - Thus we have 


a e = 


2 r2F(Tf,fc 

M 


V U 3 ~ U 1 V U 3~ 


— 7r. 


(44) 


Simplification of Eq. (1441) gives 
a e = 4, 


2 r e „ 


E(Tf, k e ) E(T|,fe e ) 


y/(r* - 3M)M L y/u% - u\ y/u%-u\i 


~ 7T. 


(45) 
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For EKBH, elliptic integral parameters can be defined as 

*5 = = si n-> 

2 y r% - 6M - 2^/(r e 0 - 3 M)M 

Modulus k e has range 0 < \k e \ 2 < 1, where 


(46) 


k e = 


1 6 M -r% + 2y/{r% - 3 M)M 


4y (r® - 3M)M 
Thus, the exact bending angle for EKBH lensing is given by 


(47) 


d e = 2, 


2r n 


V(rt - 3M)M 


K(fc e )-F(T e ,fc e ) 


- 7T, 


(48) 


where F(^,k e ) = K(k e ) defines the complete elliptical integral and i ? ('F e ,fc e ) is an incomplete 
elliptical integral. 



0 5 10 , 15 20 25 

b 


FIG. 8: For extreme Kiselev black hole lensing, the bending angle a e is a function of the impact parameter 
b (setting M = 1). In this case, the bending angle also depends on the value of the quintessence parameter 
a. In this figure, a = 0.125 is the value for the extreme case while a = 0 is for the Schwarzschild black hole 
bending angle taken as a reference. 

Figure El shows that by increasing the value of b, the bending angle decreases. The dashed 
curve shows the bending angle for EKBH, while the solid curve shows the bending angle for the 
Schwarzschild black hole. Both curves display the same behavior since they have one horizon. In 
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EKBH lensing, the event horizon is twice the Schwarzschild’s horizon (rfj). However, the difference 
between these two bending angles is that in the extreme case, the bending angle is larger than 
the Schwarzschild black hole bending angle because if we increase the value of the quintessence 
parameter <r, the bending angle will also increase. 


VIII. ALTERNATIVE APPROACH FOR FINDING BENDING ANGLE 


Gravitational lensing phenomena involves the study of the null geodesic equations. When the 
solution of the space-time geometry ((2|) extends, an event horizons exist at r + and r_, see Eq. [71 
Our main interest is in the region that lies between the horizons, which is called the photon sphere 


' ps 


[Eq. 115] . Therefore, the deflection will occur when a ray of light passes through that region 


with the closest approach r a . In order to compute the bending angle a we need to compute the 
value of the impact parameter b. If we divide Eq. m with (I12j) we obtain 

1 


d<f> 

dr 




(49) 


Now, for the closest approach r = r a and ^ \r=r 0 = 0, we will have 

b(r 0 ) = 


1 -™-(JT 0 

I o 


(50) 


By substituting Eq. (l50jl in Eq. (l49l) . we obtain 


dcj) 

dr 


(51) 


) ( 1 - - cjr n ) - I 1 - ^ - ar 


2 M 


We adopt the procedure of [3j], thus we will use the following bending angle formula: 


9 r 

a = Z —dr — 7r. 
.h dr 


(52) 


By using Eq. the deflection angle for a light ray becomes 

i 

dr 


a(r Q ) = 2 


l-^L-aro) - 1-^-crr 


— 7T. 


(53) 


The geometry of a lensing phenomenon is shown in Fig^ [9j This figure is commonly called “the 


lens diagram”. The lens equation can be expressed as 


m 
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FIG. 9: The lens diagram. The positions of observer (O), source (S), lens (L) and image (I) are shown. The 
observer-lens, observer-source, and lens-source distances are represented by Dol, Dos and Dls, respectively. 


tan (5 = tan 6 


£>ls 

Ads 


tan(a 


6) + tan 0 


(54) 


where .Dls is the distance from the lens to the source and Dos is the distance from the observer 
to the source. We also have 


b(r 0 ) = -Dql sin#, 


(55) 


where Dol is the distance from the observer to the lens. Angular positions of source and images 
are represented by (3 and 6, respectively while the deflection angle due to a black hole is denoted 
by a as it is shown in Figure EB 

Now, if we convert the distance and the impact parameter in terms of the Schwarzschild black hole 
radius, we find 




r 

2M ’ 


X 0 = — , 
2M 


b(r Q ) = 2Mb(X 0 ), 


dol — 


-Dql 


_ Dqs 
2 M 


d\ s 


Dls 
2M ' 


2 M ’ 


(56) 
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From here, we will introduce a new quintessence parameter a£ = 2Mo in terms of the Schwarzschild 
radius. Using Eqs. (f55l) and ([56]) in Eqs. (1531) . (1501) . ([7]) . and (fT5l) respectively, we get 


a{X Q ) = 2 


dX 


~ 7I\ 


X„ X \l (Xo ) l 1 Xo a £ X o) (l X a t X 

X, 


b(X 0 ) = = d Q \ sin 6 


-^H = X— — \/-T — cr €, 

2ai Gt v 4 


7T ps = 


1 - VI - 3ct £ 




(57) 

(58) 

(59) 


where Xu denotes the distance from the horizons and Xp S is the distance from the photon sphere. 
In order to find the position of images, we need to solve Eq. (1541) for the source position /3 along 
with Eqs. (1571) and (l58l) . 

n 

Generally, for a circular symmetric lens, the magnification is given by [1(| 


d = 


sin / 3 d/3 


-i 


sin# dO 

Here, the tangential magnifications and the radial magnifications are respectively defined as 


Mt 


By differentiating both sides of Eq. 


sin/3^ - 1 
sin# 


fjj V 


d6 


we get 


37] 


d£ = / coSj8 \2r _ db r / 

d6 V cos 6) - d OR \ V 


cos 9 \ 2 / da 


— - i)) 

do ) I 


mos(a — 0)J \d0 

where ^ By taking the derivative of Eq. ([571) with respect to X 0 , we obtain 


(60) 


(61) 


(62) 


da 

dXo 


X 2X 0 -3-G e X* 


/*. 2X o 


dX. 


(63) 


(x/) “ X~ 0 - a e X °) - - X - a i x ) 

Finally, by differentiating Eq. (1511) with respect to 6 on both sides and doing some simplifications 
we get 


dX 0 X o 1 - £ - OlXo) 2 x/l - (^) 2 (1 - i - °lXo) 


i-1 


d9 


2d a \ ( 2X ° 3 


oi x i) 


(64) 


IX. WEAK FIELD LIMIT 


We are going to take some approximations in this section. If the source and the lens are 
aligned, then we can approximate tan (3 ~ f3 and tan# ~ #. For the relativistic images we can write 
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A a = 2nir + A a n , (where n is an integer) and 0 < Aa„ < 1. Hence, we can replace tan(a — 9) by 
Aa n — 9. If the ray of light reaches the observer after it turns around the black hole, the deflection 
angle a must be very close to 27 t. Therefore, Eq. (154(1 becomes 


P = 9- =L Aa„ = 9-^- A a n , 

-L'OS ^os 


(65) 


and the impact parameter is b = d Q \9. 

Relativistic images are formed only if the ray of light passes very close to the photon sphere. 
For the closest approach X 0 , it is convenient to write 


X 0 — X ps + e 


0 <£< 1 . 


For a Schwarzschild black hole, the approximated deflection angle will be [11] 

-2 + V3 


a ~ —2 In 
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£ — IT. 


Therefore, we shall also look for a similar approximation [if 


( 66 ) 


(67) 


a = — Ahx(Be) — n, (68) 

where A and B are positive numbers that we take from [stJ. However, in our case these numbers 
will depend only on ap. Therefore, we will have 


lim [-(x„-A'p «)*'“* 


X 0 ^A' p 


dx n 


B = lim 

A 0 —»A ps 


exp | ^ — aexa J^ t+7r ^ | 


X 0 - Xp s 


(69) 


Now, by taking the value of X Q from Eq. (1661) and by putting that expression into Eq. ([58]). we 
get the impact parameter in terms of £ as 

X ps + £ 


b(X ps + £) — 


1 - 


-Xps+£ 


— a A X ps + £ 


(70) 


If we use a Taylor expansion in £ up to second order in Eq. (1701) . we get the impact parameter as 


b = C-D£ 2 


(71) 


where 


C = 


D = 


1 - VI - 3 a 1 


<J i\/ — 1 T \/l — + 2at 

<Ji | - (2 - 2 V 1 - 3 a^j + < 7 ^ 8 - 5V1 - 3(7^ - 6<r^j 
2^ — 1 + VI — + 2a \J—2 + 2V1 — 3<r 1 + 5 an — 2 aey/T— 3ai 


(72) 

(73) 
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Now we can find the value of e from Eq. (ED using 


£ = 


C-b 


D 


(74) 


Finally, If we substitute the value of e [Eq. (1771) into (f68l) ]. we obtain the approximated bending 
angle expression 


a = — Ain 


B 


c - d ol e 


D 


— 7T. 


(75) 


X. RELATIVISTIC IMAGES 


Virbhadra and Ellis defined “relativistic images” of a gravitational lens as those images which 
occur due to light deflections by angles a > 3 tt/2 [10[. Similarly, when (3 = 0 and d > 2-7T, the 
location of relativistic “Einstein rings” are specified [36j]. For a fixed value of (3 , we can get 0 
related to the positions of corresponding images. Thus, we can do an approximation using a first 
order Taylor expansion around a = 2ni: for the position of the nth relativistic image 


9 ~9°- Pn&Otn, 


(76) 


where 9 = 6° at a = 2nir and 


- 

Pn — ) | a=2mr 

da 


(77) 


For the value of 9 we take Eq. ED, and we get 

1 


6 = w 
d Q \ 

no = J_ 
n doi 


c ~ ^ ex p{^(« + 7r )} 


C — 


B 2 

D 


_| exp {_i( 2 n + i M 


Taking derivatives in (1791) and then substituting into (1771) , we obtain 


(78) 

(79) 


1 r 2D 


Pn — 


d Q i Lab 2 


exp 


{t( 2 " + 1 M 


(80) 


From Eq. ()76|) . we have 




On ~ 9° 

Pn 


Using Eqs. (1791) and (f80l) in (l8ll) . we get 
A r - r2 




r doiB' z ( 2 
exp 


2 LI D 




(81) 


(82) 
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Substituting Eq. (18T1) into (l65l) yields 


13 = 9 n - ^-Aa r 

^OS 


(83) 


Putting Eq. 

P = 


into ([83]) . we get 



(AB 2 

r 2 / \ i 11 



r B 2 C r 2 ( \ i -ii 

1+ d 

L U-OS 

l 2D 6XP 

ii( 2n + 1 M}] 

Un i 

^os 


D eXP {x4( 2 "- + HI 


In order to obtain the approximate position for the relativistic images, we neglect the number 1 
because ( ^J rfo1 1) in this approximation. Therefore, we have 


On = 


d n 


2D 


“ p { ~i( 2n + x )'}] 13+ i l c - w exp (it ( 2n+ 0 


d ls d Dl iab 2 ^i at 1 v" j r 1 d Dl r b 2 ^i a r" 1 v'■ jj• (85) 

Here in Eq. (1851) . if the source, lens, and image are perfectly aligned then j3 = 0 and we can obtain 
the Einstein ring with angular radius 


9 " = ^[ c -# exp 


{t( 2p+1 M 


The amplification of the nth relativistic image is given by 

/3 d/3 -i 

T n W n 

Tangential magnification for relativistic images is 


d 0 


_ w n _ ^os 

= 4doi L AB 2 


2D 


exp 




+ 


f3d 0 \ 


c-A exp 




Radial magnification for relativistic images is 


d@ri 


d n 


Hr — 


2D 


■ exp 




( 86 ) 


(87) 


( 88 ) 


(89) 


d/3 d\ s d Q \ IAB 2 

Thus, the total amplification of the nth relativistic images can be calculated by combining both 
tangential magnification Eq. (1881) and radial magnification Eq. (1891) in (1871) . which yields 


Hn = 


1 d 0 s 

|/3| di s doi 


2D 

AB 2 


exp 


{t( 2 " +1 M 


- r 1 

- Ldoi l 


^ C ~ ^2 eX P 




(90) 


Here, if the observer, lens and source are aligned (/3 = 0), the amplification will diverge. Therefore, 
the size of the relativistic images will become very small and the brightness will be low. For the 
total magnification of relativistic images, the sum of the relativistic image is taken into account 


,, _ nyoo _ ‘3dos 00 

HR n=lHn — |0| ^ ^ J n=l”nrn- 


(91) 


Now, by using the geometric series E£‘L 1 a n = for |a| < 1, the total magnification of the 
relativistic images will be 


HR 


2 d os 2D 
W\ disd 2 os AB 2 


D 

B 2 l 


/■ 


exp (—12 k/A\ 

r/ A ) 


1 — exp — 8ir 


M 


exp I — 67 t 


1 — exp —4-7T 


,/a) 


(92) 
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XI. DISCUSSION 

We have studied the GL scenario for nonextreme, naked singularity and extreme cases for 
KBH. We discussed the null geodesics for these three cases in order to study the behavior of the 
scalar held. We observed that effective potential and the null-geodesics trajectories depend on the 
quintessence parameter. From Figs. Q] and 0 we found that the potential does not have a minimum 
value so there is no stable circular orbit for photons. Moreover there are only unstable orbits for 
all cases. We also studied the behavior of the light in the lensing process of KBH. From Figures [2] 
and [71 we ensured that as the value of impact parameter b is increased the value of ro increases. 
We have worked with the quintessence held, so due to the effect of quintessence parameter a, the 
situation gets reversed i.e., closest approach ro decreases by increasing the value of b and light goes 
closer to the KBH. Moreover, when a reaches to 0.125, the ro remains constant with respect to b. 
For this, we calculated the equation of the path and the bending angle a. After that, we converted 
this expression in terms of elliptic integrals. The bending angle depends on the value of a. For 
each case, a has different limits. We solved the elliptical integrals numerically and studied their 
behavior via plots in Figs. [3 HI El and 0 

We also studied a GL phenomenon for nonextreme KBH (0 < a < %jj)- In this case, it can be 
seen from Fig. 0 that as the value of the impact parameter increases, the bending angle decreases. 
Nevertheless, for the whole process, for large value of a, light goes closer to the black hole and the 
bending angle would be larger. Furthermore, when we compared it with the Schwarzschild case, 
we observed that a s is smaller than the bending angle for the nonextreme case. 

For a GL phenomenon for EKBH, we have a = From Fig.0 we noticed that as the impact 
parameter b increases, the bending angle d e for EKBH decreases. When we compared it with 
Schwarzschild black hole, we observed that its behavior is similar to the Schwarzschild black hole 
bending angle a s and nonextreme bending angles, since EKBH has only one horizon which is twice 
the Schwarzschild’s horizon. However, a e is greater than the a s . 

To study GL phenomena for naked singularity, we took a > gjjj. In this case, the behavior of 
the light is totally different as there is no horizon and the value of the closest approach r Q will 
remain constant with respect to b. From Fig. HI it can be seen that as we increase the value of 
6, the bending angle increases. However, from Figs. HI and 0 one can conclude that the bending 
angle is smaller for large a. For the case of naked singularity, we found that the bending angle is 
larger than the nonextreme, extreme and Schwarzschild cases. (The order of the bending angles is 
naked singularity > extreme KBH > nonextreme KBH > Schwarzschild black hole) Additionally, 
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the behavior of a naked singularity bending angle is almost opposite both nonextreme KBH and 
extreme KBH bending angles. We calculated the bending angle by another approach in Sec. 
VIII and we found that the results are similar for both approaches. We have also calculated the 
approximated bending angle by using the weak field limit. The expression for the magnification of 
relativistic images are also derived. 

One can generalize this analysis and comparison for the Reissner-Nordstrom black hole sur¬ 
rounded by quintessence matter and the study of relativistic images can also be done more rigor¬ 
ously. This type of work might be important for studying the highly redshifted galaxies, quasars, 
super massive black holes, exoplanets, dark matter candidates and so on. 
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